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ABSTRACT 

The Atacama Cosmology Telescope (ACT) is currently observing the cosmic microwave background with 
arcminute resolution at 148 GHz, 218 GHz, and 277 GHz. In this paper, we present ACT's first results. Data 
have been analyzed using a maximum-likelihood map-making method which uses B- splines to model and 
remove the atmospheric signal. It has been used to make high-precision beam maps from which we determine 
the experiment's window functions. This beam information directly impacts all subsequent analyses of the 
data. We also used the method to map a sample of galaxy clusters via the Sunyaev-Zel'dovich (SZ) effect, 
and show eight clusters previously detected with X-ray or SZ observations. We provide integrated Compton-}? 
measurements for each cluster. Of particular interest is our detection of the z = 0.44 component of Abell 3128 
and our current non-detection of the low-redshift part, providing strong evidence that the further cluster is more 
massive as suggested by X-ray measurements. This is a compelling example of the redshift-independent mass 
selection of the SZ effect. 

Subject headings: cosmic microwave background - cosmology: observations - galaxies: clusters: general - 
methods: data analysis 
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1. INTRODUCTION 

A new generation of experiments is measuring the cos- 
mic microwave background (CMB) at arcminute resolu- 
tions. Within the p ast year alon e, results from the 
South Pole Tel e scope ([Staniszewski"e t al. 2Q QJ), ACBAR 
(iReichardt etal 2QQ9a). AMiBA fometsu et alJ l2009l) . 
APEX- SZ (Reichardt et al. 2QQ9b) . the Cosmic Background 
Imager ([Sieyers et al.ll2009|) , t he Sunyaev-Zerdovi ch Array 
(ISharp et al.l2009l) . and QUaD (iFriedman et al.l2009l) have re- 
vealed the ^arcminute structure of the CMB with higher pre- 
cision than ever. The angular power spectrum of temperature 
fluctuations at these scales (£ > 1000) will further constrain 
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models of the early universe. Furthermore, secondary features 
such as the Sunyaev-Zel'dovich (SZ) effect and gravitational 
lensing probe the growth of structure. 

With its first science release, the Atacama Cosmology 
Telescope (ACT) now adds to these endeavors. A 6- 
meter, off-axis Gregorian telescope, it was commissioned 
on Cerro Toco in northern Chile in October of 2007. Its 
current receiver is the Millimeter Bolometer Array Cam- 
era (MB AC), containing three 32x32 arrays of transition 
edge sensor (TES) bolometers observing at central frequen- 
cies of 148 GHz, 218 GHz, and 277 GHz, with beam full- 
widths at half-maxima (FWHM) of 1.37^ 1.0^, and 0.9^, 
respectively. It has operated for two seasons and is cur- 
rently in its third season. In 2007 one month of science 
observations were made using only the 148 GHz array. The 
other two frequencies were added for the 2008 season, which 
lasted about 3.5 months. The tel escope optic a l desi gn is 
described in iFowler et al.1 (l2007h . iHincks et all (l2008h and 
ISwitzer et al.l (l2008l) report on the telescope performance 
and provide an overview of hardware and software systems. 
The MB AC design and detai ls of T ES detector propertie s 
and readout are in Niemac d (120061). iMarriage et all (120061). 
iBattistelH et"aD (I2OO81) . iNiemack et al.l (l2008l). ISwetz et al.1 
(l2008h . lThornton et al.1 (l2008h . and lZhao et al.l (l2008h .^ 

ACT is located at one of the premier sites for millimeter 
astronomy because of the high altitude (5200 m) and the dry 
atmosphere. The precipitable water vapor (PWV) had a me- 
dian value of 0.56 mm during the nights of our 2008 season. 
Nevertheless, atmospheric emission remains the largest sig- 
nal external to the receiver in our raw data, a reality for any 
ground-based millimeter-wave telescope. The atmospheric 
power dominates only at low temporal frequencies and this 
is the main reason we observe while scanning our telescope 
in azimuth. Though much of the atmospheric power is below 
the frequency of our 0.0978 Hz scans, on typical nights the 
atmosphere dominates the detector noise up to about 2 Hz. 

In this paper we present a map-making method designed to 
model and remove the atmospheric signal in a manner which 
is unbiased with respect to the celestial signal. The method 
currently produces its best results on small scales (< 1°), so 
it is well-suited to making maps of objects with small angular 
sizes. One of the most useful applications has been the study 
of our instrumental beam with high signal-to-noise maps of 
planets. The beam profile affects all aspects of data analysis, 
including calibration, and we provide the beam characteristics 
in this paper. Additionally, we present new SZ measurements 
of eight known clusters. 

We proceed as follows: in ^we introduce the map-making 
method, showing both the theory and some qualitative prop- 
erties; O describes how we analyzed our beams, and presents 
the key measured parameters along with beam maps and ra- 
dial profiles; window functions are derived in ^ shows 
a selection of clusters imaged with the mapper; and we con- 
clude in ^ 



for interpolation, discussed more below. The technique com- 
putes maximum-likelihood estimates of both the celestial and 
the atmospheric signals, using all available detectors in a sin- 
gle frequency band. We refer to it hereafter as the Cottingham 
Method. 

In the following subsections, we give a mathematical de- 
scription of the Cottingham Method ( ^2. lb . followed by a dis- 
cussion of its benefits and a comparison to the "destriping" 
method developed for Planck, which has close similarities 
( S2.2I) . Our approach for including th e eff ects of spatial vari- 
ability across the detector arrays is in ^2.31 We discuss the use 
of B- splines in ^2.41 and finish by outlining our implementa- 
tion of the method ( ^2.5\ and map-making steps (^ ^2.6t . 

2.1. The Algorithm 

The measured timestream d is modelled as a celestial signal 
plus an atmospheric component: 



d = m+Bcx + n, 



(1) 



where the pointing matrix P projects the celestial map m into 
the timestream, 5 is a matrix of basis functions with ampli- 
tudes a which model the temporal variation of atmospheric 
power, and n is the noise. The timestream of measurements 
d may be a concatenation of multiple detectors if they have 
been properly treated for relative gain differences. Through- 
out this paper, this is the case: all working detectors from one 
frequency band are processed simultaneously. 

We seek a and in, estimates of the atmospheric amplitudes 
and the celestial map, respectively. Eq. [T] prescribes that we 
subtract the atmospheric term to obtain the map estimate: d' = 
d-BoL. The maximum likelihood estim ator is then giv en by 
the standard map-making equation (^e. g . . iTegmarkll 1 997h : 

ih = {P^N-^Py^ P^N-^d' = n {d-sa) , (2) 
where we call 



n = (p^N-^p) p^N 



(3) 



the projection matrix and N = (nn^) is the noise covariance. 
The projection matrix 11 is designed in such a way that the 
map estimate is not biased, in the sense that the error, m-m, 
does not depend on m. 

Given a set of basis functions B, the Cottingham method 
minimizes the variance of the map pixel residuals with respect 
to the amplitudes a. The residuals are the differences between 
the celestial signals measured in the timestream and the map 
estimate: 

Ad = d'-Pm=d-Ba-PU {d-Ba) 

= (l-PU){d-B5c), (4) 

where 1 is the identity matrix. We differentiate x^- 



2. THE COTTINGHAM MAPPING METHOD 

In this section, we present a techniq ue for remov i ng the 
atmospheric power first described by 'Cottingham' (' 1987h 
and used by iMever et al.] (1199 Ih . Houghn et al. ( 1992), and 
iGanga et al.l ([l993^ The temporal variations in atmospheric 
signals are modelled using B-splines, a class of functions ideal 

^ Reprints of all the references in this paragraph may freely be downloaded 
from: http://www.physics.princeton.edu/act/papers.html 



= -2B^ (1 -PUfN-^ (l-PU) {d-Ba) 

= -2B^N-^ (1 -PU) {d-B5c) . (5) 

The last equality can be obtained by expanding II to its con- 
stituent elements (c.f. Eq.O and simplifying. If we define the 
following: 
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H = B^N-^ (1 -/^n) , = = Ed, 



(6) 



then when we set the derivative in Eq.[5]to zero, we have the 
simple expression: 



05 = 0. 



(7) 



This is a linear equation which is straightforward to solve for 
the atmospheric basis function amplitude estimates a. These 
can then be used in Eq. [2]to estimate the map. In fact, both 
OL and in are the maximum likelihood estimators of the atmo- 
sphere and celestial map, respectively, for a given set of basis 
functions B. We show this explicitly in Appendix lAl 

There is an arbitrary overall offset to the computed Bol 
which must be estimated to remove the background from 
maps. We return to this point in 



2.2. Discussion 

The chief strength of the Cottingham Method is that it es- 
timates the atmospheric power in a way that is unbiased with 
respect to the map estimate itself. This important but subtle 
point is encapsulated in the term (1 -PJl) in Eq. [51 whose 
effect is to project out the map estimate from the data. There- 
fore, the solution to Eq.|7]is not sensitive to the estimated ce- 
lestial temperature, but only to a time- varying term which is 
represented by the atmospheric estimate Bol. This should be 
contrasted with high-pass filtering or fitting a slowly-varying 
function to the timestream to remove low-frequency power. 
Such approaches require masking of high signal-to-noise ce- 
lestial objects (like planets or clusters) and/or multiple iter- 
ations to prevent corruption of the maps. They also remove 
low-frequency power without regard to its origin, the major- 
ity of it atmospheric, but also inevitably containing celestial 
signal, and therefore require extensive simulations to under- 
stand the effects of the filters on the final maps. 

The Cottingham Method has close similarities to the "de- 
stripi ng" technique deve l oped in particular for Planck anal- 
vsis (iDelabrouilld Il998l: iBurigana et aP 1 19991: iMaino et al.l 
l2002h . In fact, the linear algebra pres ented in g2.1lis identi - 
cal to some versions of destriping (e.g. jKeihanen et a020Q4l) . 
The destriping techniques are intended primarily to remove 
1 / f in strumental noise — thus, for example, Keihanen et al. 
(l2QQ5h i mpose a prior on th e estimate Bex based on detec- 
tor noise. I Sutton etaP (l2009l) also consider the effects of im- 
posing a prior on the atmospheric noise. On the other hand, 
we use the Cottingham Method to remove atmospheric power 
with a flat prior. A distinct feature of our method is that 
we process multiple detectors simultaneously since the at- 
mospheric signal is common across detectors (see, however, 
^2.31) . Further, our approach differs in that it u ses B- splines 
as the basis for modelling the atmosphere ( ^2.41) . 

2.3. Spatial Structure in the Atmosphere 

The Cottingham Method as presented thus far assumes that 
the atmospheric signal Bol is common among all the detec- 
tors. In fact, we know that there is also spatial structure in the 
atmosphere, meaning that in principle, each detector might 
see a different atmospheric signal. In practice, the finite tele- 
scope beam sets a lower limit on the spatial scale. We find that 
the atmospheric signal is coherent across a quarter to a third of 
the array, or about 5-7 ' . For reference, our 148 GHz channel, 
which has a 1.37' FWHM in the far-field ( ^3. 2b . is sensitive 
to an angular size of approximately 10' at a 1km distance. 



roughly the distance to a typical turbul ence layer in the atmo- 
spher e when pointed at 50° in altitude (iPerez-Beaupuits et aP 
[2{)05h . 

To account for this, we divide the 32 x 32 detector array into 
9 square sub-arrays of roughly equal size and fit for nine sep- 
arate temporal atmospheric signals BsOLs, with the subscript s 
denoting the sub-array. These can all be done simultaneously 
if we adapt Eq. [TJ 



d=Pm+S 



/Bi 
B2 
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= Pm + SB'oL' +n, 





/ai\ 
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\Ot9 / 



(8) 



where 5" is a book-keeping matrix that remembers from which 
sub-array each measurement in d came. The Cottingham 
Method proceeds exactly as before, except that we change 
B^SB' anda ^ ol' . 

2.4. The B-Spline as a Model of Atmospheric Signal 

We follow ICottinghamI (Il987h in choosing cubic basis B- 
splines for the basis functions B. B- splines are widely used 
in the field of geometrical modeling, and numerous text- 
books cover them (e.g., Bojanovetal. 1993; de Boor 2 0011: 
ISchumaked2007l) : here we summarize basic properties. Basis 
B-splines are a basis of functions whose linear combination is 
called a B- spline. The basis B-splines are fully determined by 
a knot spacing Tk and a polynomial order p\ a B- spline is flex- 
ible on scales larger than Tk, while on smaller scales it is rela- 
tively rigid. The basis B-splines bj^p(t) of order p are readily 
evaluated using the Cox-de Boor recursion on the polynomial 
order. For m knots {tj} with j = to m - 1 : 



^j,o(0 = 



if tj <t < tj+i 
otherwise 



'-'j-b,^_,(t)^-'^''''' 



^ -bjn^p-i{t\ (9) 

with j values restricted so that j-\-p-\-l <m-l. For m knot 
times, m+p-l basis B-splines cover the interval between the 
first and last knot time. The individual basis B-splines bj^p(t) 
are compact functions, such that the B- spline receives support 
from no more than p of its bases at any point. For model- 
ing the atmospheric signal, we always choose knots uniformly 
spaced in time and use p = 3 (cubic). 

Due to their flexibility on large scales, B-splines are ideal 
for modelling the slowly varying atmospheric signal. The fre- 
quency fk below which power will be removed is determined 
by the knot spacing Tk. Empirically, we find: 

A^l/2r,. (10) 

Figure \T\ shows an example of atmospheric estimation us- 
ing the Cottingham Method. The B-spline knot-spacing is 
Tk = 0.25 s, chosen for this example because it has fk = 2 Hz, 
the approximate frequency at which the atmospheric power 
meets the detector noise level. Longer knot spacings produce 
qualitatively similar results, except that they cut off at lower 
frequencies, as per Eq.fTOl 

The Cottingham Method is effective at suppressing atmo- 
spheric contamination, but some covariance between atmo- 
spheric and celestial map estimates remains. This is typically 
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Fig. 1. — An example of the Cottingham Method. The fit is done using 
300 s of data from 605 148 GHz detectors. The PWV was 0.8 mm, about 
0.25 mm higher than the median in 2008. The knot spacing is 0.25 s and 
the order is cubic. In both plots, the original signal is plotted with a solid, 
light line, the B -spline atmosphere model with a dashed line and the signal 
minus the model with a solid, dark line. Each plot has been smoothed with 
a 5-sample boxcar filter for readability. The temperature units are with re- 
spect to a Ray leigh- Jeans spectrum. Top: A portion of one of the detectors' 
timestreams. Bottom: The spectral densities for the same single detector. A 
Welch window was applied before computing the Fourier transform. 

at harmonics of the scan frequency 0.1 Hz), as exempli- 
fied in the bottom panel of Fig. [T] For this work we have 
used the white noise approximation for the detectors (N =1), 
which results in maps of bright sources which are clean down 
to the -40 dB level in most cases (©• The small resid- 
ual atmospheric-celestial covariance is manifested as striping 
along lines of constant altitude since with our 7° peak-to- 
peak, 1.5° s~^ azimuthal scans (or 4.47° at 0.958 °s~^ when 
projected on the sky at our observing altitude of 50.3°), the 
knot spacing (r^^ = 1 .0 s for beam maps (© and 0.5 s for clus- 
ter maps (©) corresponds to an angular scale smaller than 
the scan width. When noted, we fit straight lines to rows of 
pixels in the map, after masking out any bright source, and 
subtract them. We call this process "stripe removal". In both 
our beam analysis and our cluster studies, we have done tests 
which show that the bias introduced by this process is not 
significant — see ^3.21 and ^5.31 Nevertheless, future exten- 
sions of the Cottingham Method would benefit from the full 
treatment of the noise covariance. 

2.5. Implementation 

Before making maps with the Cottingham Method, some 
preprocessing must be done. The data, which are sampled 
at 400 Hz, are divided into fifteen-minute time-ordered data 
(TOD) files and the preprocessing is performed on each in- 
dividual TOD — a future paper will describe the steps which 
we only summarize here. The data acquisition electronics' 
digital anti-aliasing filter as well as measured detector time 



constants are deconvolved from the raw data. Low frequency 
signal due to cryogenic temperature drifts are measured with 
dark detectors (i.e., detectors uncoupled to sky signal) and re- 
moved from signal detectors; a sine wave with period 10.23 s 
is also fit to each timestream and removed to reduce scan- 
synchronous contamination. Calibration to units of power 
uses nightly load curves obtained by sweeping through de- 
tector bias voltages and measuring the response. Relative 
gain imprecisions are removed by using the large atmospheri c 
signal itself to flat-field the detectors (e.g.. iKuo et al.l 12004]) : 
t his i s done independently in each of the nine sub-arrays (see 
^2.31) . Finally, calibration to temperature units uses measure- 
ments of Uranus, which we estimate to have a net error of 
6%. (The beam maps require no calibration to temperature — 
in fact, the temperature calibration is obtained from them.) 
The timestreams require no further preprocessing. 

To improve the speed of the Cottingham algorithm, we ex- 
ploit the fact that the map pixelization used for calculating 
the atmospheric signal (Eqs. [5]-[7]) need not be the same as 
the map-making pixelization. In general we only use a se- 
lection of the possible pixels on the map; additionally, we 
down- sample the number of hits in each pixel. We call the 
former "pixel down- sampling" and denote the fraction of re- 
tained pixels np\ the latter we term "hit down- sampling" and 
denote the fraction of retained hits nu. Consequently, the frac- 
tion of total available data used is npXnu. Each of these 
down-samplings is done in an even manner such that there 
are no large gaps in the remaining timestream. 

We have specified four parameters for the Cottingham 
Method: the knot- spacing r^^, the pixel size ^, the pixel down- 
sampling fraction, np, and the hit down-sampling fraction, nu. 
Of these, we always choose ^ = 18'' (about 1 /3 the 277 GHz 
beam size). 
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Fig. 2. — Top panel: The ratio L of low-frequency power to white noise 
after removing the atmospheric signal, as defined by Eq.[TT] All data points 
are from the same TOD using all working detectors; only the parameters Tk, 
Up and Uh (see text) were varied to obtain each point. The x-axis is a product 
of the resamplings Up and Uh, and the colors are different knot spacings T]^ as 
indicated in the plot key. Bottom panel: The computation time required for 
the data shown in the top panel. (The segmented lines are an artifact of how 
the data were recorded.) As the fraction UpUh of total data used increases, 
the efficacy L of power removal flattens out and adding more data does not 
significantly improve the fit but only takes more computation time. In this 
example, there were ^ 103000 possible map pixels with an average of 287 
hits per pixel. 

To evaluate the effect of varying the other three variables, 
we define a figure of merit which compares the average de- 
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tector spectral density below IHz to the white noise level, 
calculated in the range 5-25 Hz: 



TABLE 1 

Summary of Beam Parameters 



1 



p25Hz , 



(11) 



where the sum runs over the A/j detectors used for the Cotting- 
ham calculation and G/ is the spectral density of the detec- 
tor after removing the estimated atmosphere signal. Fig. [21 
shows a plot of measured values of L for a selection of knot- 
spacings and pixel down-samplings. As expected, shorter 
knot spacings remove more power: note however that only 
the Ty^ = 0.5 s and r^^ = 0.25 s are capable of removing power 
up to the 1 Hz for which L is defined (c.f Eq.fTOb. Because the 
curves flatten out as UpUh increases, at a certain point adding 
more data does not substantially improve the fit. This sup- 
ports our conclusion that we only need to use a fraction of the 
data to estimate the atmosphere. 

The timing data in the lower panel of Fig. [2l were measured 
on a 64-bit Intel Xeon® 2.5 GHz processor. Computation 
time is dominated by the calculation of the variables in Eq.[6l 
In general, these go linearly with the number of data points 
d and quadratically with the number of basis functions B\ in 
the case of B- splines, the compact support of the bases can 
be exploited so that the quadratic rate is subdominant to the 
linear, as the plot shows. 

The number of basis B-splines is small enough that it is ac- 
tually feasible to solve Eq.|7]exactly. For most cases, h owever, 
the conjugate gradient method (e.g. JPress et al.lll992l p. 83ff) 
is much faster and yields indistinguishable results. Therefore, 
we use the latter in our implementation. 

2.6. Map -Making 

Once the atmospheric model Bex has been calculated for 
a TOD with the Cottingham Method, we create its celestial 
map. Maps are made in (Aa, AAcos(<2)) coordinates, where 
Aa and AA are the distances from the altitude and azimuth of 
the map center.^ We use a pixel size of 10. 6'^ per side, about 
20% the size of the 277 GHz beam; note that this is differ- 
e nt fro m the pixel size ^ used for calculating the atmosphere 

We make the white noise approximation for each detector 
and weight it by the inverse of its variance in the map estimate. 
The detector variances are obtained iteratively: we make a 
map with equal detector weights and measure the variances 
of individual detector maps against the total map, remake the 
map with the new variances and repeat until the total map 
variance converges. The atmosphere estimates returned by 
the Cottingham Method have arbitrary offsets, which can be 
different for the nine sub-arrays we use (see ^2.31) . Thus, in 
the same iterative process, we also fit for the sub-array offsets 
and remove them when coadding detectors. 

Coaddition of TOD maps is done after all of the steps de- 
scribed above. The maps are weighted by their inverse vari- 
ances (calculated after masking bright sources or clusters). 

Finally, we mention that the software used for the results 
in this paper has a completely independent pipeline from our 
main map-making software which solves for the full survey 
area coverage. It has been especially useful for studying and 
optimizing the signal extraction in small, targeted regions, and 
has provided important double-checks for our other pipeline. 

^ This is a very good approximation to the Gnomonic projection for the 
small map sizes we use. 





148 GHz 


218 GHz 


277 GHz 


Map Properties { ^3.1\ 
#TODs 
Stripe-removal? 
Map Radius () 
Mask Radius (0^ 


16 
no 
21 
18 


15 
no 
15 
9,11,13 


11 
yes 
15 

6,8,10,12 


Beam Centers ( ^3.2\l 
Major FWHM (0 
Minor FWHM (0 
Axis Angle (°) 


1.406 ±0.003 
1.344 ±0.002 
62±2 


1.006 ±0.01 
1.001 ±0.003 
137 ±9 


0.94 ±0.02 
0.88 ±0.02 
98 ±13 


Ow Wing Fits r.4l2l 
Fit Start, Oi C) 
Fit End, $2 (0^ 
Best-fit Ow (0 


7 
13 

0.526 ±0.002 


5 

7-11 
0.397 ±0.01 


4.5 
6-10 
0.46 ±0.04 


Solid Angles (£3 
SoHd Angle (nsr) 
% interpolated 


218.2±4 
2.8 


118.2±3 
4.3 


104.2 ±6 
7.2 


Beam Fits (^4l\ 
Oo (0 


0.2137 


0.1562 


0.1367 



Note. — See text for definitions of these parameters and how they are 
measured. Values for the 148 GHz and 218 GHz bands are obtained from 
the coadded, unprocessed maps whereas the 277 GHz values are average 
values from stripe-removed maps with the mask sizes indicated in the 
table. 

^ The 218 GHz and 277 GHz beam properties are averaged from the re- 
sults at these mask radii — see text.^ The fit ranges for the 218 GHz and 
277 GHz band are varied along with the mask radii so that O2 is never 
larger than the mask — see text. 

3. BEAM MAPS AND PROPERTIES 

Understanding the telescope beams, or point- spread func- 
tions (PSF's), is of primary importance for the interpretation 
of our maps since they determine the relative response of the 
instrument to different scales on the sky and are central to cal- 
ibration. For act's measurement of angular power spectra, 
the Legendre transform of our measured beam profile, called 
a window function, determines the response of the instrument 
as a function of angular scale. 

Planets are excellent sources for measuring the telescope's 
beam because they are nearly point sources and are brighter 
than almost any other celestial object. The best candidates 
for the ACT are Saturn and Mars; of the rest, Jupiter is too 
bright and saturates the detectors, Venus is available too near 
to sunrise or sunset when the telescope is thermally settling, 
and the others are too dim for exploring the far sidelobes of 
the beam. (However, Uranus is useful as a calibrator.) The 
beam maps presented in this section are from observations 
of Saturn, which was available from early November through 
December of 2008. 

3.1. Data Reduction 

Maps were made for each night-time TOD of Saturn, using 
the Cottingham Method with ry^ = ls, ^^ = 0.32 and nh = 036. 
We had more TODs than were needed to make low-noise 
beam maps, and we excluded about 1/3 of the maps which 
had higher residual background contamination, manifested by 
beam profiles that significantly diverged from those of the 
cleaner maps. The map sizes and number of TODs per fre- 
quency band are sho wn in Table [T] 

In the analysis of ^3.21 below, we sometimes compare un- 
processed maps to stripe-removed maps. Stripe-removed map 
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have been treated as outlined in ^2.41 while unprocessed maps 
have not been altered after map-making except for the sub- 
traction of a background level. The planet is masked out 
before calculating the mean map value which estimates this 
level. The mask sizes for the three arrays, whether used for 
stripe removal or background estimation, are listed in Table[T] 
In each frequency band, the selected TOD beam maps were 
coadded. Weights were determined from the rms of the mean 
background level, calculated outside the mask radius. Rel- 
ative pointing of individual detectors was measured to sub- 
arcsecond precision using the ensemble of Saturn observa- 
tions. The overall telescope pointing was determined from 
each planet observation prior to map-making and used to cen- 
ter each TOD map, so recentering of the maps was unneces- 
sary before coaddition. 



3.2. Beam Measurements 

Fig. [3] shows coadded beam maps for the three arrays using 
a color scale which highlights the features in the sidelobes. 
The 148 GHz and 218 GHz maps have striking similarities, 
most notably along the altitude (or vertical) direction where 
both exhibit more power near the top of the map. This is 
due t o the off-axis design of the telescope (see iFowler et"an 
12007*): since the 148 GHz and 218 GHz arrays sit at about the 
same vertical offset from the center of the focal plane, their 
resemblance along this axis is expected. Note that we recover 
structure in these map at a < -40 dB level. The 277 GHz map 
is clearly inferior, showing residual striping in the scan direc- 
tion, although this occurs below -20 dB. We believe this is 
from a combination of the brighter atmosphere at 277 GHz, 
as well as detector noise correlation induced by large opti- 
cal loads (such as Saturn), which we are still investigating. 
Nonetheless we are still able to measure the 277 GHz solid 
angle to about 6% (see below), and work is underway to im- 
prove it. 

For the 148 GHz and 218 GHz arrays, we do our beam anal- 
ysis on maps which have not had stripe-removal because this 
process removes the real vertical gradient from the maps (see 
Fig.O. Nonetheless, the solid angles (see below) from stripe- 
removed maps are within la of the values from unprocessed 
maps. On the other hand, the larger residual striping in the 
277 GHz maps necessitates the use of stripe-removed maps. 

The beam center is characterized by fitting an elliptical Airy 
pattern — the function describing the beam of an optical sys- 
tem with a perfect aperture — to the top ~3 dB of the beam 
map. This provides a measurement of the location of the beam 
center, its FWHM along the major and minor axes of the el- 
lipse, and its orientation, which we define as the angle of the 
major axis from the line of zero altitude relative to the beam 
center. The uncertainties in these parameters are determined 
using the bootstrap method (iPress et al.ll 19921 pp. 69 Iff) and 
give errors consistent with the standard deviation of values 
measured from the individual TOD maps. The FWHM and 
angles are listed in Table [T] They are included for reference 
but are not used in any analysis. 

We denote the beam map by B{6^ </>), where we use coordi- 
nates with radial distance from the beam center and polar 
angle (j). By definition, 5(0, The symmetrized beam is 

averaged around the polar angle: 



Another quantity of interest is the accumulated solid angle, 
which measures the total normalized power within a given ra- 
dius: 



27r 



e'de'B{e',(i)'). 



(13) 



The beam solid angle is Q.a = ^(^ = tt). 

Fig. [4] shows measured beam profiles and accumulated solid 
angles for the three arrays. We measure the beam profiles 
down to about -45 dB. If the beams exactly followed an Airy 
pattern, these data would account for 98% of the solid angle. 
Since systematic effects could corrupt our maps at the largest 
radii, we seek a way to robustly estimate the last few percent 
of the solid angle on each beam. The method is to extrapolate 
the data with a fit to the asymptotic expression for the Airy 
Pattern: 



b\0 > Of) = 



(14) 



where Of is the beam FWHM and Ow defines the wing_scale. 
EQ.fT4l is good to better than 1% beyond about 50f (ISchroeded 
l200a §10.2b).^ Knowledge of Ow allows us to infer the 
amount of unaccounted solid angle beyond the map bound- 
ary. A simple integration shows that the solid angle beyond a 
radius Oh is: 



03 

nw(0>0t) = 27r^. 



(15) 



\0). 



(12) 



We can also use this expression to estimate the amount 
of true beam power which was "mistakenly" included in the 
measurement of the background level outside the mask radius 
and subtracted from the map. In our analysis of the beam 
profiles and solid angles (including those displayed in Fig. [4]), 
we use the fits of Ow to calculate this missing power and add 
it back into the map. A new Ow is then calculated from the 
corrected map; after two such iterations the Ow fit converges. 

Fig. |4] includes over-plots of the wing estimates from the 
best-fit values of Ow on unprocessed maps. We denote the 
radii between which the fits were performed as O1-O2, and 
choose Oi ^ 50f for each array. For 148 GHz, we obtain good 
fits for any choice of 6>2 up to 13^ or about the -40 dB level 
in the profile. Thus, we use O2 = 13\ for which we fit with 
of 40 for 35 degrees of freedom. The fits to the other profiles 
are not as robust: 218 GHz has a reduced-x^ of 2.8 for O2 = 
1' and 277 GHz has reduced-x^ of 25 for 02 = 6'. Larger O2 
gave poorer fits. Consequently, for these profiles we calculate 
6>vv at different mask sizes, as indicated in Table [T] At each 
mask size we varied O2 in 2' increments, always keeping it 
lower than the mask size. The average value from the whole 
ensemble of fits gives us Oy^ and we take its standard error as 
the uncertainty. Although Eq. [141 may be too simple a model 
for these profiles, contributions to the solid angle at these radii 
are only a few percent of the total solid angle, which has an 
uncertainty dominated by the contribution of the beam at radii 
less than O2 — see below. The values of 0^; for all three beam 
profiles are listed in Table [TJ 

^ In full generality Eq. [T4]is also proportional to cos(7rZ)^/A - 37r/4), 
where D is the telescope aperture diameter and A the wavelength; we have 
smoothed over cosine cycles. 
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Fig. 3. — Beam maps for the three frequency bands: from left to right, 148 GHz, 218 GHz and 277 GHz. The maps are from coadded observations on 1 1-15 
nights (see Table [T) and have radii of 21^ (148 GHz) and 15^ (218 GHz and 277 GHz). Maps are normalized to unity and contours are in decrements of -10 dB. 
The color scale has been chosen to highlight the fact that we have made < -40 dB beam measurements of our 148 GHz and 218 GHz bands. Even in the inferior 
277 GHz map, striping is still below -20 dB. The circles show the sizes of the beam FWHM for each band (see Table A Gaussian smoothing kernel with 
(7 = 0.54^ has been applied to highlight large-scale structure; smoothing is not otherwise performed in the analysis. No stripe removal has been done on these 
maps. 
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Fig. 4. — The beam profiles (Eq.[T2) and accumulated solid angles (Eq.[T3j for the three arrays, calculated from coadded maps (see text). The beam profiles are 
shown for both unprocessed maps, with dark errorbars, and stripe-removed maps, with light errorbars. Over-plotted on each profile is the best fit of Ow (Eq.[T4j to 
the unprocessed beam profiles. The error on the profiles are standard errors from the azimuthal average. The accumulated solid angles are from the unprocessed 
maps (without any solid angle extrapolation via Eq.[T5l for 148 GHz and 218 GHz, and from the stripe-removed map for 277 GHz. Saturn is bright enough that 
the rms power from the CMB falls below all the points in these plots. 



Our Ow fits allow us to calculate precise solid angles. At 
radii smaller than O2, we integrate the normalized power in the 
map (c.f. Eq. [13]). Beyond 62, we use Eq.[T5]to extrapolate 
the remaining solid angle. (In the case of the 218 GHz and 
277 GHz solid angles, we use the smallest O2 and the largest 
mask size in the ranges shown in Table [T] Other choices from 
these ranges do not significantly alter the results.) Finally, in 
the approximation that Saturn is a solid disk, it adds half of its 
solid angle Qs to the measured instrument solid angle — this 
is shown in Appendix |Bl Thus, the total solid angle is: 



= n(0 < O2) + nw(0 > O2) - ns/2. 



(16) 



During the period of our observations, Saturn subtended solid 
angles from 5.2 to 6.0 nanosteradians (nsr). We use the mean 
value of 5.6 nsr. 



Determining the rest of the uncertainty in the solid angle is 
not straightforward since systematic errors dominate. For our 
total error, we add the estimated uncertainties of each of the 
terms on the RHS of Eq. [161 in quadrature. The uncertainty 
from Saturn's solid angle we take to be 1 nsr, both because of 
its varying angular size and to account for any systematic er- 
ror due to the disk approximation."^ The uncertainty of ftw is 
derived from the error of the fitted Ow- For ft(0 < O2), which 
dominates, we estimate the error by looking at the distribu- 
tion of values from the individual TOD maps which comprise 
the coadded map. We did this in two ways. First, we cal- 

^ The rings of Saturn add a layer of complication to its solid angle cal- 
culation, particularly since they have a different temperature than the disc. 
The ring inclination was low during our observations (< 6°) and we have 
estimated that their contribution is negligible within the error budget. 
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culated the mean and standard deviation of the solid angles 
measured in each individual map. This also reassures us that 
the coaddition step does not introduce any systematic error 
through, for example, pointing misalignments or changes in 
telescope focus from night to night. Second, we used the 
bootstrap method to generate 1000 coadded maps with ran- 
dom subsets of individual maps and used this ensemble to es- 
timate the 68th percentile (i.e., 1 a) of solid angles. These two 
error estimates were consistent with each other. 

The solid angles and their uncertainties are reported in Ta- 
ble [T] The formal uncertainties have been doubled and we 
quote them as la, in case there are systematic effects for 
which we have not accounted. In particular, the maps used for 
power spectrum estimation will come from an independent 
pipeline and will treat the instrumental response in slightly 
different ways — for example, by weighting detectors differ- 
ently. We expect the beam uncertainties to decrease as our 
analysis evolves. 

4. WINDOW FUNCTIONS 

The statistics of the CMB are frequently characterized by 
an angular power spectrum Q : 

AT(h) = ^^aifnYimin); {a},^,aim) = ^^'^^m'mQ, (17) 

where AT(h) is the CMB temperature at position h and Y£m is 
a spherical harmonic. In spherical harmonic space, the beam 
is encoded in a window function W£ describing the response 
of the experiment to different multipoles £, such that the total 
variance of a noiseless power spectrum is: 



Var = 2^^ — Q 



47r 



(18) 



In the case of a symmetric beam, the window function is the 
square of t he Legendre transform of the beam radial profile 
(IWhite & S rednicki 1995; Bond 1996): 



^ / b\0)Pmd(cose). 



(19) 



4.1. Basis Functions 

For calculation of the window function and its covariance 
we model each beam with a set of basis functions which is 
complete but not necessarily orthogonal: 



b'(0) = Y,anK(Oy 



(20) 



n=0 



Because the bea m is truncated by a cold Lyot stop 
(iFowler et al.l l2007l) . its Fourier transform is compact on a 
disk, which suggests that a natural basis with which to de- 
compose the Fourier transform of the beam image is the set 
of Zernik e polynomials that fo rm an orthonormal basis on the 
unit disk (iBorn & Woljll999h . The Zernike polynomials, ex- 
pressed in polar coordinates p and (p on the aperture plane, 
are: 

C(P,<P)=c(py""", (21) 

where m and n are integers such that n>0^n> \m\ and ^- \m\ 
is even. In the case of an azimuthally symmetric beam, we 
need only consider the m = radial polynomials, which can 
be expressed in terms of Legendre polynomials, Pn(x) as: 



(22) 



The radial Zernike polynomials have a convenient analytic 
form for their Fourier transform: 



R2n(0) = j I 



pdpe-'P'RUp) = {-irJ2nMO)/0, (23) 



where 7„ is a Bessel function of the first kind. Motivated by 
this, we adopt 



(24) 



'"'''-[To) 



as our set of basis functions to fit the radial beam profile.^ 
Here, we have introduced a fitting parameter, 6>o, to control 
the scale of the basis functions. 

4.2. Fitting Basis Functions to the Beam Profile 

Below Oi (c.f. Table[T]), we fit the bases bn of Eq.[24lto the 
measured beam profile, and beyond 6i, we use the power law 
defined in Eq.fT4lwith the parameters 6>vv listed in Table[Tl We 
assume vanishing covariance between the power law and the 
basis functions as they are fitted to independent sets of data 
points. 

We employ a nonlinear, least- squares method to solve for 
the coefficients a„ and their covariance matrix C^^ . The algo- 
rithm uses a singular value decomposition to determine if the 
basis functions accurately chara cterize the data an d also com- 
putes a goodness-of-fit statistic (iPress et al.lll992L §15.4). As 
inputs to the fitting procedure we are required to specify the 
scale parameter, 6>o, and polynomial order, ^max- We searched 
the {6^0, ^max} parameter space until a reasonable fit was ob- 
tained that kept fZmax as small as possible. For all three bands, 
^max = 13 gives a reduced ^ 1- No singular values is found 
for any of the fits. The parameters we use for each fre- 
quency band are listed in Table [T] 

4.3. Window Functions and their Covariance s 
Given the amplitudes a„ of the radial beam profile fitted to 
the basis functions and the covariance matrix C^^ between the 
amplitudes am and a„, the beam Legendre transform is: 



bi = ^anbin' 



(25) 



n=0 



and the covariance matrix of the beam Legendre transforms 
bi and bi' is: 



(26) 



and therefore the covariance for the window function is: 



(27) 



In Fig. \5\ we show the window functions for each of the 
three frequency bands with diagonal error bars taken from 

^ It may be asked why the Airy pattern, which was somewhat suitable 
for the high-^ fit in ^3.21 is not used here. We find that at low 0, it is a 
poor fit since the optics are more complicated than the perfect-aperture model 
assumed by the Airy pattern. 
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Fig. 5. — Normalized window functions (top) and diagonal errors (bottom) 
computed from the basis functions for each of the three frequency bands. 
The window functions have been normalized to unity at = 0. In practice, the 
normalization will take place over the range of multipoles corresponding to 
the best calibration. Only statistical errors are shown. 

the covariance matrix, E^^,. We observe that the window 
function for each of the frequency bands has fallen to less 
than 15% of its maximum value at ^ = 10000. The statis- 
tical diagonal errors are at the 1.5%, 1.5%, and 6% levels 
for the 148 GHz, 218 GHz, and 277 GHz bands respectively, 
as shown i n Fig. [51 They are computed following Eq. 17 of 
iPage et all (l2003h . The off-diagonal terms in the beam co- 
variance matrix are comparable in magnitude to the diagonal 
terms. Singular value decompositions of E^^, yield only a 
handful of modes with singular values larger than 10~^ of the 
maximum values: 5 modes for 148 GHz, 4 for 218 GHz, and 
2 for 277 GHz. Thus, the window function covariances can 
be expressed in a compact form which will be convenient for 
power spectrum analyses. 

For 277 GHz we estimate a 10% systematic uncertainty 
from destriping. Another source of systematic error in the 
window functions arises from the beams not being perfectly 
symmetric. The symmetrized beam window function gen- 
erally underestimate s the power in the beam (e.g.. Fig. 23, 
iHinshaw et al.ll2007h . I n pr actice, the scans in our survey field 
are cross-linked (see S5.ll) so that the quantity of relevance 
to the power spectrum calculation is the cross-linked window 
function i.e.. 



1 
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w, = — / b{£)b\t)dc^,, 
27r Jo 



(28) 



where is the polar angle in spherical harmonic space and 
h{t) is the transform of the beam profile. To study the mag- 
nitude of this effect on the window function we computed 
the fractional difference between the window function derived 
from the Legendre transform of the symmetrized beam and 
the cross-linked window function, derived for a typical cross- 
linking angle of 60°. The difference between the two was 
found to be at the 1%, level for the 148 GHz and 218 GHz 
arrays, and at the 4% level for the 277 GHz array. 

5. SZ GALAXY CLUSTERS 

In addition to beam maps, the Cottingham Method mapper 
has been used for making maps of SZ clusters. The maps and 
analysis presented in this section are the first results from ACT 
on SZ science. For this first overview, we present results from 
only the 148 GHz band, the most sensitive during our 2008 
season. We focus on known clusters, including four new SZ 
measurements. 



5.1. Data 

Table [21 lists the clusters studied in this paper, including in- 
formation on t h e m aps and a summary of the results of our 
analysis ( ^5.21 ^5.31) . Fig. [6l shows the cluster maps and com- 
panion difference maps (see below). 

Apart from planets, ACT has done no targeted observations 
of specific objects, so the cluster maps come from our reg- 
ular survey data, which were taken at two different central 
azimuth pointings, one on the rising sky and the other on the 
setting sky. Therefore, the maps presented here are "cross- 
linked", i.e., they consist of data taken with two distinct an- 
gles between the azimuthal scan direction and the hour angle 
axis. The integration time is short, ranging from about 3 to 1 1 
minutes — see Table [2l 

The clusters were discovered in a full- survey 148 GHz map 
produced by our main map-maker. (Full maps will be shown 
and discussed in a future publication.) A Wiener filter was 
constr ucted using the polytropic model of iKomatsu & Selja^ 
(1200 ih as an SZ template, and included detector noise, CMB 
power and point source contributions in the noise model. 
Clusters were then identified from the filtered maps. We make 
two points about the clusters presented in this paper: first, al- 
though our template-based detection method has some built- 
in bias, the detections presented here are significant (> 3cr in 
the filtered survey map); and second, we have only included a 
sample of our significant detections. 

Cluster map s are made using the procedures outlined in 
^2.5l and ^2.61 The knot spacing was r^^ = 0.5 s and the down- 
sampling fractions were Up ^ 0.42 and rih ^ 0.40. All m aps 
are 0.4° in diameter. Straight-line stripe removal ( ^2.6\i has 
been performed, using a 6' radius mask over the cluster decre- 
ment. 

We have made companion "difference" maps for each clus- 
ter from the same data. For each of the rising and setting 
observations, a map made from the first half of the nights' 
data is subtracted from the second half. The rising and setting 
difference maps are then coadded to produce the full, cross- 
linked coadded difference map — the same procedure used for 
the signal maps. 

The map noise, listed in Table[2l is the rms of the map com- 
puted outside a 6' mask and converted to an effective pixel 
size of one square arcminute. By examining the power spec- 
tra of the maps we found that the rms values we quote are 
dominated by the white noise level and do not have signifi- 
cant contributions from residual low-frequency power. 



5.2. Characterizing the Detections 

The cluster center positions are determined by finding the 
coldest point in the map smoothed with a 2' FWHM Gaus- 
sian kernel. The one exception is ACT-CL J0509-5345 (SPT- 
CL 0509-5342), which has a complex structure. In this case, 
we choose a center which gives a maximal signal-to-noise 
(see below). As a rough guide, we also quote the cluster 
depths, ATsz, from these smoothed maps in Tabled but we 
stress that these values should not be used for quantitative 
analysis. All other cluster properties are measured from un- 
smoothed maps. 

We quantify the significance of the detections by a signal- 
to-noise (SNR) measurement. Denoting the signal map ntxy 
and the difference map dxy, we define the SNR within an aper- 
ture size as: 



1 The ACT Collaboration 



TABLE 2 

Selection of SZ Clusters Detected by ACT 



ACT Descriptor Catalog Name J2000 Coordinates^ rms^ ^int'^ SNR(ef ATsz^ 10^^x7(6)^ 

[I^K] [min] [^K] 



RA Dec. e<2' 6> < 4' 6* < 6' 

(±0.2) (±0.6) (±1.2) 



ACT-CL J0245- 


-5301 


Abell S0295 


02^45^^28' 


-53°01'36'' 


44 


10.1 


15.2 (6.80 


-250 


0.89 


2.36 


3.91 


ACT-CL J0330- 


-5228 


Abell3128 (NE) 


03^30^^50' 


-52°28'38'' 


49 


10.3 


12.8 (4.30 


-260 


0.94 


2.69 


4.34 


ACT-CL J0509- 


-5345 


SPT-CL 0509-5342 


05^09^^20' 


-53°45'00'' 


47 


10.1 


7.7 (5.20 


-70 


0.33 


1.07 


1.50 


ACT-CL J0516- 


-5432 


Abell S0520 


05hl6"^3P 


-54°32M2'' 


55 


6.8 


4.2 (4.10 


-110 


0.19 


-0.11 


-0.55 


ACT-CL J0546- 


-5346 


SPT-CL 0547-5345 


05^46^^35' 


-53°46'04'' 


46 


9.5 


13.9 (5.80 


-250 


0.91 


2.36 


3.67 


ACT-CL J0638- 


-5358 


Abell S0592 


06^38^^46' 


-53°58M0'' 


55 


7.5 


8.1 (3.10 


-230 


0.70 


1.40 


2.07 


ACT-CL J0645- 


-5413 


Abell 3404 


06^45^^29' 


-54°13'52'' 


59 


9.3 


2.8 (2.00 


-120 


0.12 


-0.18 


-0.69 


ACT-CL J0658- 


-5556 


IE 0657-56 (Bullet) 


06^58^^33' 


-55°56M9'' 


80 


3.4 


12.1 (2.70 


-510 


1.60 


2.95 


3.56 



^ Position of the deepest point in 2^ FWHM Gaussian smoothed map, except for ACT-CL J0509-5345 which has a position which gives a maximal 
SNR (see text).^ Map rms measured outside a 6' mask and reported for a one square arcminute area.'^ Integration time, defined as the approximate 
total time (in minutes) that the telescope was pointed in the map region.^ Maximum signal-to-noise ratio (Eq. [29l and the radius at which it was 
obtained.^ Cluster depth, as measured in a 2' FWHM Gaussian smoothed map at the listed coordinates; intended as a guide to the magnitude of the 
decrement.^ See Eq.[32]and following discussion. 



SNR(0) = min < 



{mjcy±djcy) 



(29) 



where the sum is over pixels which are within a radius d from 
the cluster center, ISlid) is the number of pixels in the sum, and 
and crj are the variances of the signal and difference maps, 
respectively, calculated outside a 6' mask. Because the sign 
of a difference map is arbitrary (we subtracted one half of the 
data from the other half), we take the smaller of the possible 
two SNR. The purpose for including the subtraction is that 
coincident flux in the difference map indicates that the signal 
map probably contains spurious flux. 

In Table [21 we quote SNR(6i) at the value of 6i which max- 
imized the SNR. 

5.3. Integrated Compton-Y Values 

The SZ effect occurs when CMB photons inverse 
Compton- scatter off hot ele ctrons in clusters of galaxies 
(ISunvaev & ZeldovichlflQTQh . The imprint on the CMB is 
proportional to the integrated electron gas pressure: 



AT 

TcMB 



kB(JT 



nieC^ 



(30) 



where the integral is along the line of sight, m^. He, and Te 
are the electron mass, number density, and temperature, re- 
spectively, (Jt is the Thomson scattering cross-section, and 
the variable y is the Compton-j parameter. The function f{x) 
encodes the dependence on frequency: 

f{x)= [xcoth(x/2)-4] [l + fe(x,r,)], (31) 

with X = hv/kBTcMB- For hot cl usters, a relativi stic correc- 
tion fe(x, r^) may be required (iRephaelil Il995h . in which 
case / must be included in the integral of Eq. [SOlin the non- 
isothermal case. However, in this paper we assume non- 
relativistic / for simplicity. 

A robust measure of the SZ signal is the integrated 
Compton-j parameter, since it is model-indepen dent and sim- 
ply sums pixels in the maps (iBenson et al.ll20Q4l) : 



Y{e) = jj d^e'y{0'\ 

\e'\<e 



(32) 



where 6 is the angular distance from the cluster center. We 
use steradians as the unit of solid angle, so Y is dimension- 
less. As an example, it is plotted for ACT-CL J0638-5358 
in the lower panel of Fig. [71 the upper panels show its az- 
imuthally symmetrized temperature profile. The values of Y 
at 2^4^ and 6' are shown for each cluster in Table [2 The er- 
rors were determined by calculating the standard deviation of 
Y values from maps of random patches of the CMB, made in 
the same way as the cluster maps. We find la uncertainties of 
0.2x10-1^ 0.6x10-10 and 1.2x10-^^ for 7 at 2^ 4^ and 6^ re- 
spectively. These values are dominated by systematic errors, 
including uncertainties in the background level, contributions 
from CMB in the map, and residual contamination from the 
stripe removal. Accordingly, they should not be interpreted as 
the significance of the detection; the SNR values (see Table[2l) 
serve that purpose. 

Note that in two clusters— ACT-CL 105 16-5432 
(Abefl S0520) and ACT-CL J0645-5413 (Abefl 3404)— we 
measure negative values of YiA') and Yi6'). Their maps show 
that the measured SZ signal is compact and the negative 
values are consistent with noise. 

As a check that the choice of knot spacing (r^^ = 0.5) is 
not creating a significant bias via covariance of the celes- 
t ial s ignal with the low-frequency atmospheric estimate (see 
^2.41) . we created maps with Tk from 0.155-10 1.5^ for ACT- 
CL J0245-5301 and ACT-CL J0638-5358. The temperature 
profiles for the latter are plotted in the middle panel of Fig. [71 
Even the shortest spacing does not produce a profile which 
is significantly different from the others. Its knot spacing, 
Tyj: = 0.15 s, is the only one from the ensemble with a corre- 
sponding angular scale smaller than the map size. We con- 
clude that the results are not biased by having knots of too 
high a frequency. 

5.4. Comparisons to Previous Measurements 

The clusters shown in this paper are previously known X- 
ray, optical, and/or SZ clusters; all are massive systems. For 
four of the sources (Abefl S0295, Abefl 3128 (NE), Abefl 
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Fig. 6. — Cluster maps made using the Cottingham Method at 148 GHz, paired with their difference maps (see ^5.U . The coordinates are J2000 right ascension 
(hours) and declination (degrees). The color bars are }xK (CMB); note that the scale is different for each cluster. The gray disc in the top corner of the the signal 
plots is 2.43^ in diameter, the FWHM size of the beam convolved with the Gaussian smoothing kernel which was applied to these images. In each difference plot, 
a cross shows the coordinates of the darkest spot in its corresponding signal map (except for SPT-CL 0509-5342 (ACT-CL J0509-5345) — see text). 
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Fig. 7. — The radial profile (top/middle) and integrated Compton Y(0) 
values (bottom) for the SZ decrement of ACT-CL J0638-5358 (Abell S0592). 
The profile data are averages from the maps in 22^^ -wide annuli, and Y(0) is 
the sum of the pixe ls within a radius 6, converted to the unities s Compton- j 
parameter (Eqs.[3Qland[32l. The top panel shows the profile of the signal map 
and difference map. The middle panel compares profiles for maps made with 
different knot spacings Tk, showing that only for very short spacings is the 
profile noticeably different from the r = 0.5 s profile used for cluster analysis, 
and then not significantly so. In all of the panels, the profile centers were 
determined by the minimum of the map after smoothing with a 2' FWHM 
Gaussian profile. (The profiles were calculated from the unsmoothed map.) 

S0592, and Abell 3404), these are the first reported SZ detec- 
tions. In this section, we briefly review measurements from 
the Hterature to provide context, and point out some of the 
contributions that our new measurements make to this body 
of knowledge. 

Relevant parameters from the literature are listed in Table[3j 
references for these values are included below. Typical errors 
on Lx are small (< 20%), while those on the inferred mass are 
more substantial (~50%). Temperatures are measured values 
from X-ray spectra. We use a flat ACDM cosmology with 
Qm = 0-3 and Hq = 70kms~^Mpc~^ Masses are quoted in 
units of M500, defined as the mass within a radius having a 
mean mass density (p) 500 times greater than the critical den- 
sity, i.e., (p) = 500 X 3if 7(87rG). In the following we briefly 
discuss the clusters in the order in which they appear in Ta- 
ble [2I with the exception of the South Pole Telescope (SPT) 
clusters which are discussed together near the end of this sec- 
tion. 

5.4.1. Abell S0295 



AbeU S0295 first appeared in lAbeU et all (Il989h in their ta- 
ble of supplementary southern clusters (i.e., clusters that were 
not rich enough or were too distant to satisfy the criteria for 
inclusion in the rich nearby cluster catalog). It was also found 
to be a significant X-ray source in th e ROSAT All Sky Survey 
(RASS) Bright Source Catalog ( Vog es et al.lll999l) The spec- 
trosco pic redshift of Abell S0295 was obtained by lEdge et al.l 
(11994)) . who also reported the discovery of a giant strong- 
lensing arc near the brightest cluster galaxy. Efforts to de- 
tect the SZ effect at 1.2-mm and 2-mm with t he SE ST were 
attempted, un successfully, by Andreani et al.l (Il996h . ASCA 



observations (iFukazawa et al.l 



luminosity. The cluster mass M500 was then estimated from 
the luminosit y-mass (specifically LvfO. 1-2 .4 keV) vs. M500) 
relations from lReiprich & Bohringed (l2002h . 

5.4.2. Abell 3128 (NE) 

Until quite recently the north-east (NE) component of Abell 
3128 was believed to be part of the Horologium-Reticulum 
supercluster at z = 0.06. The X-ray morphology is clearly dou- 
ble peake d with the two com ponents separated on the sky by 
some 12 MRose etall(l2002h estimated the virial masses of the 
two components assuming the redshift of the supercluster and 
obtained a value for each of ~1.5 x 10^^ Mq. Fig. [8] shows 
our SZ measurement with overlaid X-ray contours. 



Abell 3128 




03h31m 03h30m 
Right Ascension - J2000 

Fig. 8.— ACT-CL J0330-5228 (Abell 3128 (NE)) with overlaid contours of 
X-ray emission in black. The SZ detection is associated with the NE feature 
of Abell 3128, and confirms that it is due to a more massive, higher redshift 
cluster than that at the SW lobe — a compelling example of the redshift in- 
dependent mass selection of the SZ effect. The X-ray data come from two 
separate XMM-Newton observations (Obs Ids 0400130101 and 0400130201) 
with a total exposure time of 104ks. The two observations were mosaicked 
into a single image over the 0.2-2.0 keV. Contour values are from 1.25x10"^ 
to 1.25x10"^ photons/cm^/s/arcsec^. 



yielded values (see Ta- 



ble [3]) for average temperature and soft band X-ray flux(0.1- 
2.4 keV), from which we determined the corresponding X-ray 



Recently [Werner et"an (l2007h carried out a detailed study of 
this cluster using XMM-Newton data, which revealed a more 
distant and more massive cluster superposed on the northeast- 
ern component of Abell 3128. A significant portion of the 
X-ray emission comes from this background cluster. The val- 
ues we quote in the table for redshift. X-ray luminosity, gas 
temperature, a nd M500 corre s pond to the background cluster 
and come from I Werner et al.l (l2007h . 

The large SZ decrement (SNR = 12.8) seen in the ACT 
maps is clearly associated with the NE component where the 
^ = 0.44 cluster is. We do not detect a significant decre- 
ment from th e south western component which lies at z = 0.06. 
IWerner et al.l (l2007h estimate the temperature of the higher 
redshift cluster to be 5.14±0.15 keV, which is significantly 
hotter than that of the foreground cluster (kT = 3.36 ± 0.04 
keV). This system, therefore, is a compelling illustration of 
the mass selec tion, approximately independent of redshift, of 
the SZ effect. I Werner et al.l (|2007|) note that the temperature, 
luminosity and mass estimates of the z = 0.44 background 
cluster are all subject to large systematic errors, as the cluster 
properties depend upon the assumed properties of the fore- 
ground system. A joint X-ray /SZ/optical analysis should be 
able to better constrain the characteristics of both systems and 
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TABLES 

Summary of Cluster Properties from X-Ray and Optical Studies 



ACT Descriptor 



Catalog Name Redshift 



Da Lx(0. 1-2.4 keV) M500 kT lO^o x F2500' 

[Mpc] [10"^"^ erg s-i] [lO^^Mo] [keV] 



ACT-CL 
ACT-CL 
ACT-CL 
ACT-CL 
ACT-CL 
ACT-CL 
ACT-CL 
ACT-CL 



J0245-5301 
J0330-5228 
J0509-5345 
J05 16-5432 
J0546-5346 
J0638-5358 
J0645-5413 
J0658-5556 



Abell S0295 
Abell3128 (NE) 
SPT 0509-5342 
Abell S0520 
SPT 0547-5345 
Abell S0592 
Abell 3404 
IE 0657-56 



0.3006 
0.44 
0.36 (P) 
0.294 
0.88 (P) 
0.2216 
0.167 
0.296 



920 
1172 
1037 
906 
1596 
737 
589 
910 



8.3 
3.9 
2.2 
3.5 
4.7 

10.6 
8.2 

20.5 



0.8 
0.3 
0.4 
0.6 
0.6 
1.0 
0.7 
1.4 



6.7 ±0.7 
5.1±0.2 

7.5 ±0.3 

8.0 ±0.4 

7.6 ±0.3 
10.6±0.1 



53+0-35 

'^•^ -0.06 
Q 72+0.47 

1 01+0.86 
^•■^^-0.52 
1 87+1.23 



Note. — See 35.4l for citations to the literature from which these values were obtained. The marker (P) in the redshift column 
indicates a photometric redshift measurement. 

^ Predicted value of Y within /?2500 from the Y-kT scaling relation of lBonamente et all ( I2008D . Errors come from the uncertainty 
on the scaling relation parameters. Although we do not ha ve R25QQ values for our clusters, the Y(2') measurements listed in 
Table [2] should be roughly comparable to these — see ^5.4.81 



thereby contribute to assessing the mass threshold of the ACT 
cluster survey. 

5.4.3. Abell S0520 

This optically-ric h cluster (I Abell et al.lll989|) lies at a red- 
shift of z = 0.294 (lGuzzoetal."'1999'). It is also an X-ray 
cluster, RXC JQ516-5 430 ( Bohringer et al. 2004), and has 
been detected by SPT ('St aniszewski et all 12008 ) as SPT-CL 
0517-5430. The X-ray temperature, soft X-ray fl ux and lumi- 
nosity are based on XMM-Newton observations (IZhang et al.l 
and the mass we quote in Table [3] comes from the X- 
ray-d erived gas density and temperature profiles (IZhang et al.l 
[20Q8h . 

5.4.4. Abell S0592 

The ga laxy cluster Abell S0592 was originally detected 
optically (lAbeU et al.|[T989l) . ROSAT detected it as a bright 
source in the A ll Sky Survey and its redshift (z = 0.2216) 
was reported in Ide Grandi et al.l (Il999l) . The cluster is also 
known by its REFLEX designation of RXC J0638.7-5358 
(iBohring er eTaP 120041) . The ROSAT flux and luminosity in 
the soft X-ray band (0.1-2.4 keV) are 7.5 x 10"^^ erg cm"^ s"^ 
and 1.1 X 10"^^ ergs s"^ The X-ray spectrum of Abell S0592 
from a Chandra observation (iHughes et al.ll2009l) yields an in- 
tegrated gas temperature of = 8.0 ± 0.4 keV. The soft X-ray 
luminosity implies a cluster mass of M500 = 10^^ M©. 

5.4.5. Abell 3404 

AbeU 3404, at z = 0.167 (Ide Grandi et al.ll 19991) . is the low- 
est redshift system in the ACT SZ-detected cluster sample 
presented here. It is REFLEX cluster RXC J0645.5-5413. 
The X-ray temperature, soft X-ray flux and lumino sity are 
based on XMM-Newton observations (IZhang et al.l 12008). 
These authors also provide the total cluster mass based on the 
X-ray-derived gas density and temperature profiles. 

5.4.6. IE 0657-56 (Bullet Cluster) 

We detect IE 0657-56, the famous "Bullet" cluster, at high 
significance with a strong central decrement and large in- 
tegrated Y. Previous detections of the mm-ba nd SZ signal 
from this cluster have been reported by ACBA R (iGomez et al.l 
[200I and APEX-SZ fHalverson eLal][2008|). 

The spectroscopi c redshift of IE 0657-56 was obtained by 
iTucker et al.l ([T998h . the X -ray flux came from the Einstein 



Observatory (iMarkevi tch et a l.ll2002h. the X -rav gas temper- 
ature from XMM-Newton (Zh ang et"ani2Q06l), a nd the cluster 
mass, M500, from a study by .Zhang et al.l (l20Q8h . 



ACT-CL 0658-5557 (IE 0657-56/Bullet) 
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Fig. 9.— ACT-CL J0658-5556 (Bullet Cluster) with overlaid contours 
of X-ray emission (black) and dark matter distribution (orange). The X-ray 
contours come from an 85ks-long Ctiandra observation (Obs Id 3184) and 
correspond to the 0.5-2.0 keV band. Contour values are 4x10"^^ to 2x10"^ 
photons/cm^ /s/arcsec^. The lensing data are from lClowe et all j2007l) with 
contours running from ^c = 0.12to0.39. 



Fig. [9] shows a zoomed-in plot of our SZ map with 
X-ray contours fro m Chandra and lensing contours from 
IClowe et aP (l2007h . As expected, the SZ decrement follows 
the X-ray contours more closely than the lensing contours, 
since the collisonless dark matter is expected to be offset from 
the collisional gas in this merging system. 

5.4.7. South Pole Telescope Clusters 

IStaniszewski et al.l (l2008h recently reported blind SZ de- 
tections of four galaxy clusters. Only one of them (SPT- 
CL 517-5 430) is a previously known cluster, Abell S0520 
(see g5Al. 

The physical properties of the SZ clusters — photometric 
redshifts, luminosities, and mass estim ates — have been re- 
ported by iMenanteau & HughesI (|2P09^ based on optical and 
X-ray data. Here we summarize some of their findings (see 
TableO. All four clusters have central elliptical galaxies asso- 
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ciated with them whose luminosities are consistent with those 
of clusters in the Sloan Digital Sky Survey. Their mass esti- 
mates from their optical and X-ray luminosities also suggests 
that these are fairly massive systems. In Table [3] we list their 
M500 estimates from their X-ray luminosities. 

We make strong detections of SPT-CL 0547-5345 
and SPT-CL0509-5342 and also see Abell S0520 (SPT- 
CL 0517-5430) with moderate SNR (c.f., Table [2]). There 
are visible similarities between the ACT maps presented here 
and those from SPT. An interesting feature is the bright spot 
to the north-west of SPT-CL 0509-5342 which figures promi- 
nently in both maps and for which we could find no radio or 
infrared catalog source. 



SPT-CL 0528-5300 Difference 




Fig. 10. — A map and corresponding difference map centered on the coor- 
dinates of SPT-CL 0528-5300, a cluster candidate detected by the South Pole 
Telescope but not yet detected with our instrument. The units of the axes are 
right ascension (hours) and declination (degrees), and the color bar is ATcmb 
(loK). The map rms noise is 45 \jlK. 

We are unable to confirm the detection of SPT-CL 
0528-5300. Fig. [Tol shows a map centered on the cluster 
coordinates. There is no measureable SNR (Eq. [29]) for a 
putative decrement centered at its coordinates. Based on 
the map noise, we report a 2a non-detection at the 90iLtK 
level . Recent weak-lensin g mass estimates for the SPT clus- 
ters (iMcInnes etal.ll2009h indicate that SPT-CL 0528-5300 
has a lower mass (M500 ~ 2 x lO^^M^r^ , scaled from their M200 
assuming a typical factor of 0.6, e.g., iReiprich & Bohringerl 
(I2002h ) than any of the other known clusters that we have de- 
tected in the SZ with ACT (see Table [3]) . Further study of the 
X-ray and optical properties of SPT-CL 0528-5300 will be 
necessary to more accurately predict its expected SZ signal. 

5.4.8. Comparison with Previous SZ measurements 

Although the large masses of the ACT-detected clusters we 
report here offer strong support for the reality of our detec- 
tions, we also compare the quoted integrated Compton-j pa- 
rameters for consistency with expections from previous SZ 
cluster studies. For this we use the Y-kT scaling relation from 
iBonamente et al.l (l2008h (using values for "all clusters" from 
their Table 2). Predicted values are given in the last column 
of Table O The Y values in the scaling relation were inte- 
grated within 7?25oo, the radius where the average cluster mass 
density is 2500 times the critical density. We do not have 
precise Rism values for our clusters, but estimates of Rism 
range from about 1' to 3 ^ so the predicted values of 7(2500) 
should, to first order, be roughly comparable to our Y(2') val- 
ues. With that proviso, the predicted and measured Y values 
agree to within Icr for all but two clusters: Abell 3404, whose 
X-ray temperature predicts a much higher Y value than we 
measure, and Abell 3128 (NE), where the cluster temperature 



predicts a much lower Y value than measured. The latter is 
a complex system which could hav e a larger mass than pre- 
viously thought. Additionally, the IBonamente et al.l (l2008h 
scaling relation was measured at 30 GHz. At 148 GHz, the 
point source contamination is different. This might explain 
why a measured Y value is lower than the prediction based on 
the 30 GHz scaling relation. A larger and better- studied sam- 
ple of ACT-detected clusters will be necessary before drawing 
conclusions about scaling relations. 

6. CONCLUSIONS 

We have described a maximum-likelihood mapping algo- 
rithm which uses B-splines to model atmospheric signal and 
to remove it from the data. 

The method has been used to make high precision (< 
-40 dB) beam maps, with solid angles in the 148 GHz, 
218 GHz, and 277 GHz bands of (218.2 ± 4) nsr, (118.2 ± 
3)nsr, and (105.2 ± 6) nsr, respectively. The beam profiles 
and window functions will be important for all subsequent 
analyses of ACT's data. 

Additionally, we have made maps showing SZ detections of 
eight previously discovered galaxy clusters. Our high-a de- 
tection of the z = 0.44 component of Abell 3128, and our cur- 
rent non-dectection of the low-redshift part, corroborates ex- 
isting evidence that the further cluster is more massive. This 
is a compelling example of the redshift-independent mass se- 
lection of the SZ effect. 

The maps presented in this paper will be made 
pubHc by Jan 1, 2010 through the LAMBDA 
site (http://lambda.gsfc.nasa.gov/) and ACT site 
(http://www.physics.princeton.edu/act/). 
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APPENDIX 

A. THE COTTINGHAM METHOD AS A MAXIMUM LIKELIHOOD ESTIMATOR 
Eq.[T]can be written in the matrix form, 



d={P B) 



which has the maximum likelihood estimator: 



N-'d = 



'(PP) 


(PB)' 




(BP) 


(BB) 





where we use the shorthand notation (XY) = X^N ^Y. The inverted matrix evaluates to: 

[(PP) - (PB)(BB)-^ (BP)] -(PP)~^ (PB) [(BB) - (BP)(PP)-^ (PB)] 

-(BB)-^ (BP) [(PP) - (PB)(BB)-^ (BP)] [(BB) - (BP)(PP)-^ (PB)] 

Since this matrix is symmetric, and (XY)^ = (YX), we can rewrite the lower-left component as: 

(BB)-\BP) [(PP)-(PB)(BB)-\BP)]~^ = ^(PP)-\PB) [(BB)-(BP)(PP)-\PB)]~^y 

= [(BB) - (BP)(PP)-^ (PB)] (BP)(PP)-^ . 

Thus, the solution for the atmosphere is: 



(Al) 



(A2) 



(A3) 



(A4) 
(A5) 



a = [(BB)-(BP)(PP)-\PB)] B^N-^d - [(BB)-(BP)(PP)-\PB)] (BP)(PP)-^P^ N'^d. (A6) 

To show that this is equivalent to the solution presented in ^2.11 we observe that the definitions in Eqs.[3]and[6]of ^2. II can be 
recast: 



(d=B^N-\l-PJl)B = (BB)-(BP)(PP)-\PB), 

cl) = B^N-\l-PU)d = [B^N-^-(BP)(PP)-^P^N-^]d. 

This reduces Eq. lA6l to: 

which is the same as Eq.|7]of ^2.11 

B. A PLANET' S SOLID ANGLE CONTRIBUTION TO THE BEAM SOLID ANGLE MEASUREMENT 



(A7) 



(A8) 



Denote the instrument response with P(n) and the power emitted by the planet with Po'^(n), where Pq is the peak power emitted 
and ^ is a normalized distribution describing its shape. The coordinate w is a two dimensional vector describing the position on 
the sky, with n = at the planet center. The measured beam map, B, is the convolution of the true beam, B, with the planet: 



B(n) = 



P(n) _ // dQ.n'B(n-n')^(n') 



P(0) JJdnn'B(-n')^(n') ' 
The measured solid angle is then (c.f. Eq.[T3]): 

JJdnn'B(n-n^)^(n^) _ [// Jl],.^(fiO] [JJ dnnB(n)] 



dQn 



// dnn,B(-n')^(n') 



JJ dn„>B(n')'^(n') 



JJ dn„,B(n')^(n'y 



(Bl) 



(B2) 



where in the second equality we brought the denominator outside the outer integral, and in the numerator we switched the order of 
and then shifted the dummy variable for the integral over B. In the last equality we recognized that the integrals in the numerator 
evaluate to the solid angles of the planet and the true instrument beam, respectively. If the planet is much smaller than the beam, 
we can expand the beam appearing the integrand of the denominator in a Taylor series: 
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B(n) = l + VB{0)n + -nH(0)n + 



(B3) 



where H is the Hessian matrix of the beam. At the beam center, being the peak, the gradient vanishes. If we assume a symmetric 
beam, then n -H n = (V^5/2)|fip, and we can write: 



1 + 



Jf2„/|nf*(n') 



(B4) 



where /af is the second raw moment of the planet shape ^. In the small planet approximation we are making, the second term in 
the brackets is small. Thus: 



(B5) 



For a disk, /i^" = Vt\/2iT, and both a Gaussian beam and an Airy pattern have V^B{0) = -47r/l^A. (For the Airy pattern, this 
is easiest to see by expanding the Bessel function in a power series and differentiating.) Thus, we have the result that Qa ^ 

nA+^^/2. 
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